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The nature of the unknown sources of ultra-high energy cosmic rays can 
\ be revealed through the detection of the GZK feature in the cosmic ray spec- 

trum. The only two experiments that have probed this energy range, AGASA 
and HiRes, have apparently conflicting results. HiRes measured a flux consis- 
tent with the GZK feature while AGASA reported a larger than expected flux of 
so-called Super-GZK particles. Here we emphasize that neither experiment has 
gathered the statistics necessary for making a definitive measurement of the GZK 
cutoff. The photo-pion production responsible for the GZK feature is stochastic 
for energies around the cutoff leading to large fluctuations of the spectrum for 
low statistics measurements. We show that the results from AGASA and HiRes 
results are within about 2a of one another by simulating 400 spectra for a range 
of input spectral indices normalized to the number of events above 10 19 eV for 
each experiment. If a 15% systematic correction in energy is applied to both 
experiments, the agreement between the experiments improves considerably and 
the best fit input spectral index becomes ~ 2.6 for both data sets. Our results 
clearly show the need for much larger experiments such as Auger and EUSO, 
that can increase the number of detected events by 2 orders of magnitude. Only 
large statistics experiments can finally prove or disprove the existence of the GZK 
feature in the cosmic ray spectrum. 
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1. Introduction 

The most extreme accelerators in the universe remain a mystery. Particles 
with energies above 10 20 eV have been observed to arrive on Earth from directions 
in space that are not clearly associated with any particular source. The mystery 
of the origin and nature of these particles should be resolved once we have a 
statistically significant measurement of the spectrum and angular distribution of 
these particles. Next generation experiments such as the Pierre Auger Project 
and the EUSO experiment should be able to finally resolve this mystery. 

Presently available data show a number of interesting hints of possible 
resolutions to this mystery. The High Resolution Fly's Eye fluorescence experi- 
ment has reported hints of a flux decrease around the Greisen-Zatsepin-Kuzmin 
(GZK) feature. A significant decrease in the flux of ultra-high energy extragalac- 
tic protons was predicted by Greisen [?] and Zatsepin and Kuzmin [?] following 
the discovery of the cosmic microwave background (CMB). Protons with energies 
above 7 x 10 19 eV lose a significant fraction of their energy through the photo-pion 
production off the CMB as they traverse tens of Mpc in the universe. If the HiRes 
result is confirmed with much better statistics the nature of the highest energy 
cosmic rays are most likely to be extragalactic protons. 

The giant airshower ground array, AG AS A, does not see the suppression 
expected at the GZK feature. They reported about 11 events above 10 20 eV 
implying a flux above the expected level due to photo-pion production. The excess 
flux seems incompatible with extragalactic protons and a number of alternatives 
have been proposed to address this discrepancy [?]. In addition, the distribution 
in the sky of arrival directions reported by AGASA shows hints of clustering in 
small scales. These may be the first hints of point sources. 

Here I discuss the significance of the apparent inconsistencies between 
the two experiments and re-emphasize the lack of sufficient statistics in both 
experiments to address the very important question of the origin of the highest 
energy particles ever observed. This presentation is based on results reported in 
the recent work of DeMarco, Blasi, and Olinto (2003) [?] (see also [?]). 

2. AGASA versus HiResI 

The two largest experiments to measure the flux of UHECRs, AGASA 
[?] and HiResI [?], reported apparently conflicting results shown in Fig. IT1 The 
figure shows that the HiResI data are systematically below AGASA data and 
that HiResI sees a suppression at ~ 10 20 eV that resembles the GZK feature 
while AGASA does not. The number of events above 10 19 eV, 10 19,6 eV and 10 20 
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Fig. 1. Flux of cosmic rays multiplied by E s versus the energy, E. Circles show the 
AGASA spectrum from [?] while squares show the HiResI spectrum from [?]. 



eV for AGASA and HiResI are shown in Table H~1 We also show in the same table 
the number of events above the same energy thresholds, when a systematic energy 
decrease of 15% is applied to the AGASA data while a systematic energy increase 
of 15% is applied to the HiResI data. These systematic energy shifts are within 
the reported uncertainties in energy measurements in both experiments. (See this 
proceedings for a discussion of energy errors in AGASA [?] and HiRes [?].) 

In order to understand the difference between AGASA and HiRes data, 
we first determine the best fit injection spectrum for each data set separately 
by running 400 simulations of extragalactic proton propagation normalizing to 
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Table 1. Number of events for AGASA and HiResI detected above the energy thresh- 
olds reported in the first column. 
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Table 2. x 2 for fits to AGASA and HiResI data above 10 18 5 eV, 10 19 eV, and 10 19 - 
eV for varying spectral index 7. In parenthesis the number of degrees of freedom. 



the data column in Table El The injection spectrum is taken to be a power law 
with index 7 between 2.3 and 2.9 with steps of 0.1. For each injection spectrum 
we calculated the x 2 indicator (averaged over 400 realizations for each injection 
spectrum). The errors used for the evaluation of the \ 2 are due to the square 
roots of the number of observed events. As reported in Table 123 the best fit 
injection spectrum can change appreciably depending on the minimum energy 
above which the fit is calculated. In these tables we give xt(N), where iV is the 
number of degrees of freedom, while the subscript, e, is the logarithm of E t h (in 
eV), the energy above which the fit has been calculated. The numbers in bold- 
face correspond to the best fit injection spectrum. These fits are dominated by 
the low energy data rather than by the poorer statistics at the higher energies. 

For energies above 10 18 ' 5 eV, the best fit spectra are E~ 2A for AGASA 
and E~ 2S for HiRes. If the data at energies above 10 19 eV are used for the fit, 
the best fit injection spectrum is E~ 2S for AGASA and between E~ 2 - 7 and E~ 2 ' 8 
for HiRes. If the fit is carried out on the highest energy data (E > 10 19 ' 6 eV), 
AGASA requires an injection spectrum between E~ 2,5 and E 1 " 2 ' 6 , while E~ 2 ' 8 or 
E~ 2 ' 9 fits the HiRes data in the same energy region. Note that the two sets of data 
uncorrected for any possible systematic errors require different injection spectra 
that change with Eth- 

In order to quantify the significance of the detection or lack of the GZK flux 
suppression, in Tables I3~1 and l4~1 we report the expected (namely simulated) num- 
ber of events above the energy threshold indicated for different injection specta 
and the discrepancy between the simulations and the data in Table El (i n paren- 
thesis) in numbers of standard deviations, a, after comparison with the data. 
From Tables EU and we can see that while HiResI is consistent with the exis- 
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Table 3. Number of events expected above Eth for different injection spectra assum- 
ing the AGASA statistics above 10 19 eV. In parenthesis are the number of standard 
deviations, a, between the expected number of events and the observed number. In 
square brackets are the combined error bar between simulations and observations, 

Cfot- 
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Table 4. Number of events expected above Eth for different injection spectra assum- 
ing the HiResI statistics above 10 19 eV from Table IT1 In parenthesis are the number 
of standard deviations, a, between the expected number of events and the observed 
number. In square brackets are the combined error bar between simulations and 
observations, a tot . 



tence of the GZK feature in the spectrum of UHECRs, AGASA fails to detect 
the same feature, but only at the ~ 2.5<r level for the best fit injection spectra. 

In Figs. ED and EU we show the uncertainties in the predicted fluxes for 
AGASA and HiRes respectively, for a injection spectrum with spectral index 
2.6. These plots show clearly the low level of significance that the detections 
above E GZ k have with low statistics. The large error bars that are generated by 
our simulations at the high energy end of the spectrum are mainly due to the 
stochastic nature of the process of photo-pion production: in some realizations 
some energy bins are populated by a few events, while in others the few particles 
in the same energy bin do not produce a pion and get to the observer unaffected. 
The large fluctuations are unavoidable with the extremely small statistics available 
with present experiments. On the other hand, the error bars at lower energies are 
minuscule, so that the two data sets (AGASA and HiResI) cannot be considered 
to be two different realizations of the same phenomenon. The discrepancies at the 
lower energy range must be due to systematic errors in at least one if not both 
experiments. 

Taking into account the (theoretical) error bars determined by the sim- 
ulations, we report in Tables ED and the standard deviations due to the low 
statistics (cr sim ). The observed number of events above the same threshold also 
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Fig. 2. Simulations of AGASA statistics with injection spectra E~ 2 ^ (upper plot) 
and E~ 2 ' 8 (lower plot). The crosses with error bars are the results of simulations, 
while the grey points are the AGASA data. 
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Fig. 3. Simulations of HiRes statistics with injection spectra E 2,6 (upper plot) and 
E~ 2 - 7 (lower plot). The crosses with error bars are the results of simulations, while 
the squares are the HiRes data. 
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has a observational error bar a b s associated with that enregy range. The com- 
bined error can be estimated by a tot = [a 2 im + a 2 bs ) 1 / 2 . We show in Tables ED 
for AG AS A and |3] for HiRes the standard deviations in the simulations (a S j m ) in 
parenthesis and the standard deviation with the combined error in units of atot i n 
square brackets. The effective discrepancy between predictions and the AGASA 
data is at the level of 2.1 — 2.5a. A definitive answer to the question of whether 
the GZK feature is there or not awaits a significant improvement in statistics at 
high energies. 

As seen in Fig. [EJ the difference between the AGASA and HiResI spectra 
is not only in the presence or absence of the GZK feature: the two spectra, 
when multiplied by E 3 , are systematically shifted by about a factor of two. This 
shift suggests that there may be a systematic errors either in the energy or the 
flux determination of at least one of the two experiments. Possible systematic 
effects are discussed in [?] and [?] for the AGASA collaboration and in [?] and 
[?] for HiResI. At the lower end of the energy range the offset may be due to the 
notoriously difficult determination of the AGASA aperture at threshold while the 
discrepancies at the highest energies is unclear at the moment. In any case, a 
systematic error of ~ 15% in the energy determination is well within the limits 
of systematic errors reported by both collaborations. 

In order to illustrate the difficulty in determining the existence of the GZK 
feature in the observed data in the presence of systematic errors, we split the 
energy gap by assuming that the energies as determined by the AGASA collab- 
oration are overestimated by 15%, while the HiRes energies are underestimated 
by the same amount. The number of events above an energy threshold is again 
reported in Table IT1 The values of x 2 obtained in this case are reported in Table 

m 

For AGASA, the best fit injection spectrum is now between E~ 2 ' 5 and 
E~ 2S above 10 19 eV and above 10 19,6 eV (the x 2 values are identical). For the 
HiRes data, the best fit injection spectrum is E~ 2 ^ for the whole set of data, 
independent of the threshold. It is interesting to note that the best fit injection 
spectrum appears much more stable after the correction of the 15% systematics 
has been carried out. Moreover, the best fit injection spectra as derived for each 
experiment independently coincides for the corrected data unlike the uncorrected 
case. This suggests that combined systematic errors in the energy determination 
at the ~ 30% level may in fact be present. 

The expected numbers of events with energy above 10 19,6 eV and above 10 20 
eV with the deviation from the data are reported in Tables W\ and 171 while HiResI 
data remain in agreement with the prediction of a GZK feature, the AGASA 
data seem to depart from such prediction but only at the level of ~ 1.8a. The 
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Table 5. x 2 f° r AGASA and HiResI in which a correction for a systematic 15% 
overestimate of the energies has been assumed for AGASA and a 15% underestimate 
of the energies has been assumed for HiResI. 
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49 ±6.9 (+0.2) [+0.1] 
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Table 6. Expected number of events above when a systematics of —15% is intro- 
duced in the energy determination of AGASA. The number of standard deviations 
between simulations and observations , a, is listed in parenthesis. In square brack- 
ets are the discrepancies calculated with a combined error bar of simulation and 
observation uncertainties, a tot- 



systematics in the energy determination further decreased the significance of the 
GZK feature, such that the AGASA and HiResI data are in fact only less than 
2a away from each other. 

We can use the same procedure illustrated above to evaluate the effect 
of the error bars in the simulated data compared to the data corrected by 15%. 
The results are reported in square brackets in Tables ED (for AGASA) and E] 
(for HiRes), showing that the effective discrepancy between expectations (with 
uncertainties due to discrete energy losses and cosmic variance) and AGASA data 
is even smaller, only at the 1.5a level. 

3. Conclusions 

The statistical significance of the UHECR spectra measured by the two 
largest experiments in operation, AGASA and HiRes, is not sufficient to deter- 
mine the nature and origin of these particles. The discrepancies between the two 
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Table 7. Expected number of events above Eth when a systematics of +15% is intro- 
duced in the energy determination of HiResI. The number of standard deviations 
between simulations and observations , a, is listed in parenthesis. In square brack- 
ets are the discrepancies calculated with a combined error bar of simulation and 
observation uncertainties, a tot- 



experiments is below 2.6 a, but AGASA sees hints of a flux above the predicted 
GZK suppression and hints of point sources in arrival direction distribution while 
HiRes sees hints of the GZK feature. If confirmed by future experiments with 
much larger statistics, the increase in flux relative the GZK prediction hinted by 
AGASA would be of great interest. This may signal the presence of a new com- 
ponent at the highest energies that compensates for the expected suppression due 
to photo-pion production, or the effect of new physics in particle interactions. 

Identifying the cause of the systematic energy and/or flux shift between 
the AGASA and the HiRes spectra is crucial for understanding the nature of 
UHECRs. This discrepancy has stimulated a number of efforts to search for the 
source of these systematic errors including the construction of hybrid detectors, 
such as Auger, that utilize both ground arrays and fluorescence detectors. If 
systematic errors of about 15% are present in both experiments the two spectra 
can be shifted into agreement with a best fit injection spectrum of i? -7 with 
7 = 2.6. 

With the low statistical significance of the possible discrepancy between 
AGASA and the expected fluxes (or the HiResI data), it is inaccurate to claim 
either the detection of the GZK feature or the extension of the UHECR spectrum 
beyond Eqzk at this point in time. A new generation of experiments is needed 
in order to finally give the proper answer to this question. In Fig. ^]we report 
the simulated spectra that should be achieved in 3 years of operation of Auger 
(upper panel) and EUSO (lower panel). The error bars reflect the fluctuations 
about the mean spectrum with an injection spectrum E~ 2 ' 6 . The energy region 
where statistical fluctuations dominate the spectrum is moved to ~ 10 20 ' 6 eV for 
Auger [?], allowing a clear identification of the GZK feature. This fluctuations 
dominated region stands beyond 10 21 eV for EUSO [?]. 
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Fig. 4. Predicted spectra and error bars for 3 years of operation of Auger (upper 
plot) and EUSO (lower plot). 
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